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Abstract 

Stabilization of passively unstable thrust-levitated vehicles can require significant control inputs. 
Although thrust vectoring is a straightforward choice for realizing these inputs, this may lead 
to difficulties discussed in the paper. This paper examines supplementing thrust vectoring with 
Variable-Speed Control Moment Gyroscopes (VSCMGs). The paper describes how to allocate 
VSCMGs and the vectored thrust mechanism for attitude stabilization in frequency domain and also 
shows trade-off between vectored thrust and VSCMGs. Using an H2 control synthesis methodology 
in LMI optimization, a feedback control law is designed for a thrust-levitated research vehicle and 
is simulated with the full nonlinear model. It is demonstrated that VSCMGs can reduce the use of 
vectored thrust variation for stabilizing the hovering platform in the presence of strong wind gusts. 

Nomenclature 

Platform parameters 

m: total mass. 

J(r 7): total inertia moment. 

IjP- inertia moment of fan rotational parts about spin axis. 

/|p inertia moment of a momentum wheel about spin axis. 

A: collection of prop spin unit vector components, 

c: first moment of inertia of the total system. 
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(7(e): Direction Cosine Matrix (DCM) from the inertia frame to the body frame 

T{uj s , 6 ),t(u) s , 0): net force and torque generated by four ducted fans 

t^ s : control torque for fan speed. 

T Ci B : control torque for a momentum wheel speed. 

cib'- spin unit vector of a momentum wheel 

Momentum wheel speed. 

VSCMG parameters 

t s , T g \ wheel speed control torque and girnbal torque 

I(,: total moment of inertia about spin axis, output axis, and girnbal axis. 

/*: wheel inertia moment about spin axis. 

ICi- total moment of inertia of the i-th VSCMG. 
r/i: girnbal angle of the i-th VSCMG 

Cq : the DCM from the body frame to the initial girnbal frame of the i-th VSCMG 
Cg, (’i'll)'- the DCM from the initial girnbal frame to the current girnbal frame of the i-th VSCMG 
(sm(??i) 0 cos{t]i))Cq i 

G{7])- ( sin ^ 0 cos ^ 2 )) C G 2 

( sin(r] 3 ) 0 cos{r] 3 ))C^ 3 

_(sm(j? 4 ) 0 cos(r/ 4 ))(7g, 4 _ 

(cos(r] i) 0 - sin(rn))C , Q i 

K ( y ( cos M 0 - sin(rj 2 ))C^ 2 

(■ cos(r] 3 ) 0 - sin(r] 3 ))C^ 3 

_(cos(r/ 4 ) 0 - sin(r] 4 ))CQ 4 _ 

1 Introduction 

Hovering vehicles levitated by vectored thrust systems frequently have their center of mass located 
above the point of application of the thrust vector, resulting in unstable or neutrally stable open- 
loop attitude dynamics. Such vehicles, and their stability characteristics, have attracted attention 
recently [1, 2, 3, 4, 5] for their suitability for a number of practical applications. A natural and 
conventional choice of control for stabilizing vehicle attitude is to actively vary the thrust vector. 
Particularly for small vehicles with low moments of inertia, however, the thrust vector correction 
necessary for stabilizing a passively unstable vehicle can be quite active, and may induce unsteady 
flow effects which are difficult to model. This difficulty is important because an unstable vehicle 
requires a high-performance control system which, in turn, requires an accurate model of the 
dynamics. This paper examines using Variable-Speed Control Moment Gyroscopes (VSCMGs) 
to, at least, partially take over the role of stabilization from thrust vectoring in these vehicles. 
These devices have the advantage that it is straightforward to accurately model their behavior. 
Additionally, by taking over a portion of the stabilization activity from the thrust vectoring system, 
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the amplitude and higher frequency content of the control from the latter can be reduced, making 
it, in turn, easier to model. 

In Ref. [1], it has been suggested to enhance directional stability of the airframe by introducing 
an angular momentum bias generated using a momentum wheel. The approach can circumvent 
the problem of control-induced unsteady aerodynamics on the closed-loop stability during attitude 
stabilization. Detailed nonlinear simulations in Ref. [1] of the NASA Flying Test Platform (NFTP) 
under a strong wind gust condition showed the potential of the approach in terms of attitude 
robustness. In Ref. [2], the feasibility of using a set of VSCMGs was demonstrated to obtain con- 
trolled performance beyond what is possible through a bias momentum stability-enhanced system. 
Using a Lyapunov stability condition, a nonlinear control law was synthesized for the NFTP system 
which incorporates vector thrusting, a bias momentum wheel and VSCMGs and was simulated in 
the presence of strong wind gusts. The improved performance of the NFTP system using VSCMGs 
was demonstrated in Ref. [2]. 

This paper describes control allocation with VSCMGs and vectored thrust machinery in attitude 
stabilization of the NFTP. Since the VSCMGs can generate internal control torques to effectively 
and very reliably redistribute angular momentum to control the platform independently of ducted 
fan and control vane aerodynamics, they are assigned to generate control torque required for attitude 
stabilization of the platform in high frequency domain. The ducted fans with control vanes are 
assigned to stabilize the platform in low frequency domain in order to mitigate unpredictable 
unsteady aerodynamic phenomena. A feedback control law is synthesized minimizing H -2 norm of 
the system augmented with weighting functions. Nonlinear simulations demonstrate the trade-off 
of control efforts between VSCMGs and vectored thrust machinery (fan speed and control vanes ). 

This paper contains the following sections. In Section 2, NASA Flying Test Platform is described 
and equations of motion are derived to calculate the internal torque generated by VSCMGs and 
the effects on motion due to girnbal angle variation. In Section 3, the detailed control synthesis 
framework is described, which involves control allocations of VSCMGs and vector thrusting in 
the frequency domain. In Section 4, nonlinear simulation results are presented and this paper is 
concluded with a brief summary in Section 5. 

2 NFTP Description and Equations of Motion 

2.1 The NASA Flying Test Platform Description 

The NFTP is a testbed for exploring combining thrust vectoring and internal angular momentum 
exchanges in unstable thrust-levitating vehicles. The system, pictured in Fig. 1, is levitated by 
thrust from four electric motors driving shrouded fans. Downstream of each fan are two thrust 
vectoring vanes for a total of eight vanes. All motors and vanes are independently controlled. The 
system is currently equipped with a flywheel, aligned with the vehicle’s body vertical axis, that 
generates an internal angular momentum bias, but the vehicle is simulated with both the flywheel 
and a set of VSCMGs. The momentum wheel speed is regulated as constant with a simple PI 
feedback control law. The detailed physical parameters of the vehicle and the angular momentum 
bias are given in Ref. [2]. 

In this paper, four VSCMGs are included in the dynamics, configured in the pyramid shape in 
Ref. [6]. Each wheel inertia moment about spin axis is 0.005 Kg m 2 and the initial wheel speed 
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of each VSCMG is an user-selected parameter that scales the authority of the CMG systems. The 
initial speed is of importance, because it comprises a nominal “set point” in the control design, 
similar to the trim values of the propulsive control effectors. The total principal moment of inertia 
for each VSCMG is roughly 

I c = diag ([/£,/£,/£]) = diag([0.0053, 0.0026, 0.0026])(Kg nr 2 ) (1) 

The output torque capacity of each VSCMG [2] is roughly 10 N-m with initial wheel speed 9000 
rpm. 


2.2 Equations of Motion 


The equations of motion of the NFTP with four VSCMGs have been introduced in Ref. [1] with 
the following states: total linear momentum p. total angular momentum h, fan angular momentum 
along spin axis h a , flying wheel angular momentum along spin axis hs z , VSCMG wheel angular 
momentum hw.i along spin axis and VSCMG angular momentum he . 2 along girnbal axis. The 
total linear and angular momentum are difficult to be measured directly; therefore the equations 
of motion are rewritten in order to re-define states which can be measurable directly by physical 
sensors. The redefined states are as follows: platform velocity v, platform angular velocity uj, fan 
motor speed uj s , momentum wheel speed fl#, VSCMG wheel speed fl and girnbal angle rate 77 . 

From Ref. [1] , the equations of motion in vectrix form are 
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The momenta are related to the velocities as follows 
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In general applications with VSCMGs such as attitude control of satellites in Refs. [7, 8, 9, 10], 
the change of the total moment inertia due to VSCMG girnbal angle variations are generally small 
enough to be ignored. In the NFTP, this assumption may be not hold. Note that the total moment 
of inertia J is a function of girnbal angle rj in this paper and is written as 

4 

J(v) = Jo + Y}°G^r)Chf IcACgM)C° G } (10) 



After algebraic manipulation of the time derivatives in Eq.(6), the equations of motion of the 
platform can be rewritten as follows: 

x = M- 1 G f (x,9,r t ,e) + M~ 1 U l (11) 

where x = [v w r) G u> s £Ib] T and U\ = [0 r g t s t^ s Tq ] t . The matrix M and function Gf 
are defined as follows: 

a = MM. I (12) 
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Note that the term V appears in Eq.(14) since the total inertia moment is a function of VSCMG 

girnbal angle. The torque r c acting on the platform in Eq. (13) is 

r c = —w x (Al^ptOg + al^^B + Kj I cV + A(r/) r I\y£l) — — f] I\y& (15) 

From Eq.(15), it is observed the gyroscopic effect due to prop rotation, the moment wheel rotation, 
girnbal angle rate change, and VSCMG wheel rotation. The last term of Eq.(15) represents the 
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internal torque generated by the VSCMGs. In order to present the trade-off effects between vector 
thrusting and VSCMGs in attitude control for the platform, it is required that the internal torque 
is explicitly formulated in the equations of motion. 

It is noticed that the matrix inversion of M € 7\l 19x19 at each time step is required to simulate 
the equations of motion of Eq.(ll). It is not computationally tractable since the computation cost 
(memory and CPU time) of inverting large dimension matrices is high. Here, using the matrix 
inversion lemma [11], the matrix inversion is calculated by using matrix multiplications and the 3 
by 3 matrix inversion at each time step, which is computationally tractable. 


Kinematic equations of motion of the platform are written as 

Xi = C T (e)v (16) 

e = 0.5[0.5(1 - e T e )/ 3X 3 + e x + ee T }ui (17) 

where e is the Modified Rodrigues Parameters (MRP) [8, 9]. The direction cosine matrix C(e) is 


C(e)=I + 


4(l-|e| 2 ) 8 

(i + M 2 ) 2 (i + N 2 ) 2 


e x e x . 


(18) 


2.3 Actuator Models 

Based on bench test of a single representative ducted fan, the electric motor is modeled as 

Imiosi = t f . {u Si ) + r^. = kiu? s . + k 2 u s + t^ s .{S P wm) (19) 

where k\ is an aero-drag coefficient due to prop rotation and k 2 is back EMF resistance coefficient 
due to armature rotation. Note that the control torque . is a function of PWM signal ( Spwm ) 
generated by the motor controller (Aveox SH 96 Motor controller). The dynamics of the electrical 
loop of the DC motor is fast enough to be ignored. To validate the model, step responses of fan 
speed with predicted step responses are plotted in Fig. 2. To model the servo motor (JRDS 8411), 
frequency sweep signals of vane angle commands are applied to the servo. Deflected angles are 
measured by potentiometers in analog signal. From the input and output data, frequency transfer 
functions are obtained by using MATLAB. The transfer functions are also calculated at different 
motor speeds and are fitted by the first order model. Hereafter, the servo model is used for 
control synthesis. 


3 Control Synthesis Framework 


In this section, a feedback control synthesis frame work for the flying platform described. The 
control objectives are to make the platform robustly stable and to minimize its attitude variations 
for hovering in the presence of wind gusts. The control design problem is formulated into minimizing 
the H 2 norm of weighted command-tracking errors and state variations with penalizing control 
actions shown in Fig. 3. An H 2 control law can be designed on a linearized model, which is 
calculated using Jacobian linearization around a given trim point. A trim condition at hovering is 
defined by setting state derivatives in Eq. (11) to zeros with additional parameter setting 

v = w = fj = f) = 0. (20) 
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Here, ( ) notation denotes trim or nominal rotation speed in the case of CMG wheels. Note that 
trim condition can be satisfied at any girnbal angle of VSCMGs. Here, the trim girnbal angle is 
set to zero for convenience. It is possible to choose a trim girnbal angle based on controllability of 
VSCMG with the given VSCMG configuration. The details are, however, out of the scope of this 
paper. 

To penalize state variations, performance weighting matrices W p shown in Fig. 3 are defined as 

W p = block diag(1000/4x4, 0. 1/4x4, IOO/4X4) (21) 

to compensate different physical units of each state. In this design frame work, the states cor- 
responding the weighting matrix W p are [H lo s /)] • The weighting matrix W e to penalize the 
command tracking error signal in Fig. 3 is set to 

W e = block diag(0. 1/3x3, 500/3x3, 2/3x3, 100/3x3)- (22) 


For hovering, the commands of X , e, V and uj are set as zero. 


In this paper, an attitude control scheme with vectored thrust machinery and VSCMGs is posed 
in which VSCMGs generate internal control torque to stabilize the platform in high frequency 
domain and vectored thrust from ducted fans with control vanes generate torque to stabilize the 
platform in low frequency domain. The reason of two control actuators allocated in two frequency 
domains is to mitigate adverse unsteady aerodynamic effects that may be provided by rapid control 
vane motions in order to reject fast dynamic disturbance in the presence of strong wind gusts. In 
Ref. [1], it has been shown that rapid motion of control vane action is required to stabilize the 
platform under adverse weather condition without using VSCMGs and bias angular momentum. 
It was also shown that vane control authority was not enough to stabilize the platform under 
strong wind gusts in Ref. [1]. In order to allocate control actions of VSCMGs and vectored thrust 
machinery on each frequency domain, the weighting function W ac t shown in Fig. 3 is defined as 


= bl “ k dia *< 17 ^TT 7 ^ i^TT 7 " 4 ' <23) 


where the entries in W ac t correspond to the inputs, read from top to bottom. The frequency re- 
sponses of the weighted function in the W ac t block are shown in Fig. 4. 

In this paper, the crossover frequency of the weighting function can be defined based on the 
mathematical model of the electric motor and the servo motor model. It, however, takes an im- 
portant role in control allocation. The crossover frequency of each frequency weighting functions 
in Eq. (23) should be defined based on an unsteady aerodynamic effect bench test. In order to 
estimate the critical frequency of unpredicted aerodynamic force and moment generated by rapid 
control vane motion, a single shrouded prop unit has been bench-tested. Note that constant weight 
in the weighting function W ac t on the VSCMG’s wheel torque command is set in frequency domain 
for simplifying a control synthesis problem. After generating an augmented open loop system shown 
by the dashed in Fig. 3 with the given weighting functions, a feedback control synthesis problem is 
formulated in H 2 LMI optimization. 

To investigate the effectiveness of VSCMGs in attitude control for the platform, feedback control 
laws are designed at different the wheel trim, or initial, speeds of the VSCMGs such as H = 1000, 
2500, 5000, and 9000 (rpm), respectively, and the wheel speed of the momentum wheel is set as 
2000 rpm. 


7 



4 


NFTP Simulation Results 


4.1 Wind Gust Disturbance Model 

To simulate wind disturbances during hovering or low airspeed operations, crosswinds are assumed 
which impinge on a vertical cylinder attached to the upper part of the platform and which gen- 
erate forces and moments via drag. The resulting aerodynamic drag on this cylindrical column 
is used to provide moment disturbances from the crosswind and turbulence impinging on the ve- 
hicle. The resulting disturbance model consists of forces and torques in the pitch and roll axis. 
The wind gust velocity is decomposed in two directions of pitch and roll axis. Each component 
of the wind gust is assumed to be independent and Gaussian with mean value (8 ft/sec) and 
standard deviation (50 ft/sec). To generate a colored gust model, the 3rd order low pass filter 

( 7t 77 ) is used. The disturbance force due to the wind gust is only assumed 

v 0.001s 3 + 0.03s 2 + 0.03s + r 

as the aerodynamic drag on the column with Cp = 0.4 at the Reynolds number (5 x 10 5 ~ 10 6 ). 
The disturbance torque is simply calculated as l a _ c _ x D where l a _ c is a position vector to the aero- 
dynamic center of the column from the platform C.G. The calculated disturbance force and torque 
acting on the platform is shown in Fig. 5 and are used for all simulations. 


4.2 Simulation Results 

In this section, simulation results based on the full nonlinear model with the feedback control 
laws designed based on each linearized model at different trim VSCMGs wheel speeds. In the full 
nonlinear model, a fan speed regulator is included and the momentum wheel speed is assumed to 
be constant as 2000 rpm for all simulations. In the following simulations, the aero-drag torque 
on the momentum wheel is ignored for convenience. The command signals for position, velocity, 
attitude angle and angular velocity of the platform are set to zero for hovering under the wind gust 
disturbance acting after 1 sec in all simulations. The attitude responses of the platform are shown 
in Fig. 6 for the case: Cl = 1000 rpm and Fig. 7 for the cases: Cl = 2500, 5000, and 9000 rpm. 
The plots in Fig. 6 and 7 imply that VSCMG actuators reduce Euler angle response effectively. 
It is observed from Figs. 6 and 7 that the attitude variations are smaller with faster wheel speed 
of VSCMGs since fast wheel speed of VSCMGs can generate large internal control torque on the 
platform. The vane deflection, fan speed, VSCMG girnbal angle responses, and VSCMG wheel 
speed responses are shown in Figs. 8 through 12, to compare control efforts for each case. From 
the plots in Figs. 8 and 9, it is observed that the control effort of fan speed and vane deflection can 
be reduced by increasing VSCMG control efforts. Note, from Figs. 9 and 10, that larger values of 
wheel speed Cl in the VSCMGs reduced the girnbal angle variations and thrust vector angles needed 
for stability. Notice in Figs. 11 and 12 that the VSCMG Cl variations are less than 1.6 % of the 
initial wheel speed for all cases. It is also observed that the fast wheel speed (9000 rpm) requires 
smallest variations of the wheel speed during simulation for holding attitude. 

From the simulation results, it is also observed that the RMS (Root Mean Square) Euler angle 
responses decrease as VSCMG Cl values increase. To describe the trade-off of control efforts between 
vectored thrust (fan speed and vane deflection angle) and VSCMGs, the RMS time responses girnbal 
torque, girnbal angle rate, fan speed, wheel speed and vane deflection are shown in Figs. 13 and 
14. It is observed from Figs. 13 and 14 that control action of vector thrust machinery (RMS of fan 



speed and vane deflection angles) decreased as action of VSCMGs (RMS girnbal torque command) 
increase. Also, notice that fast wheel speed VSCMGs require small girnbal angle rate variation at 
constant girnbal torque. It corresponds to the fact that the internal control torque acting to the 
platform is proportional to the product of wheel speed and girnbal angle rate. 


5 Conclusion 


This paper has described attitude control using VSCMGs and vectored thrust machinery in the 
presence of strong wind gusts demonstrated on an experimental testbed, the NFTP. A feedback 
control law for the NFTP has been designed by minimizing H 2 norm with weighting functions which 
penalize the use of vectored thrust in high frequency domain to mitigate unsteady aerodynamic 
phenomena generated by rapid control vane motions. The required control torque for attitude 
stabilization in high frequency domain can be generated by the VSCMGs. Simulation results show 
the trade-off control effects between VSCMGs and vectored thrust for the NFTP with different 
VSCMG wheel speed levels. For the practical use of VSCMGs on aero- vehicles, the performance 
requirement under given environment is required to size VSCMGs. 


Acknowledgments 


This research of the first author’s work was supported by National Aeronautics and Space Admin- 
istration under NASA Contract No. NAS1-02117. The authors thank their colleagues Elvin Ahl 
and Jeffrey Massie for constructing the NFTP frame, and David Cox and Mario Smith for setting 
the bench test of the single ducted fan. 


References 

[1] Lim, K.B., Shin, J.-Y., Moerder, D.D., and Cooper, E.G., “A New Approach to Attitude 
Stability and Control for Low Airspeed Vehicles,” in AIAA Guidance, Navigation, and Control 
Conference and Exhibit, (Providence, R.I.), August, 16-19 2004. AIAA Paper 2004-5008. 

[2] Lim, K.B., Shin, J-Y., and Moerder, D.D., “Variable Speed CMG Control of A Dual-Spin 
Stabilized Unconventional VTOL Air Vehicle,” in AIAA 3rd Unmanned Unlimited Technical 
Conference, Workshop and Exhibit, (Chicago, IL), September, 20-23 2004. AIAA 2004-6537. 

[3] Army Aerial and Joint Services VTOL Rotocraft Goverment/Industry / Academia Workshop, 
(Williamsburg, VA), January, 27-29 2004. 

[4] Fleming, J., Jones, T., Gelhausen, P., and Enns, D., “Improving Control System Effectiveness 
for Ducted Fan VTOL UAVs Operating in Crosswinds,” in AIAA 2nd Unmanned Unlimited 
Systems, Technologies, and Operations, (San Diego, CA), September, 15-18 2003. AIAA 2003- 
6514. 


9 



[5] Lim, K.B., Shin, J.-Y., Cooper, E.G., Moerder, D.D., Khong, T.H., and Smith, M.F., “An 
Overview of The NASA Flying Test Platform Research,” in AIAA Guidance, Navigation, and 
Control Conference and Exhibit , (Austin, TX), August, 11-14 2003. AIAA 2003-5775. 

[6] Yoon, H. and Tsiotras, P., “Singularity Analysis of Variable-Speed Control Moment Gyros,” 
Journal of Guidance, Control, and Dynamics , Vol. 27, No. 3, 2004, pp. 374-386. 

[7] Hughes, P.C., Spacecraft Attitude Dynamics , ch. Cp.3. New York: John Wiley and Sons, 1986. 

[8] Tsiotras, P., “Stabilization and Optimality Results for the Attitude Control Problem,” Journal 
of Guidance, Control, and Dynamics, Vol. 19, No. 4, 1996, pp. 772-779. 

[9] Schaub, H., Vadali, S.R., and Junkins, J.L., “Feedback Control Law for Variable Speed Control 
Moment Gyros,” Journal of the Astronautical Sciences , Vol. 46, No. 3, 1998, pp. 307-328. 

[10] Wie, B., Spacecraft Vehicle Dynamics and Control. AIAA Education Series, Reston, VA: AIAA 
Inc., 1998. 

[11] Zhou, K., Doyle, J., and Glover, K., Robust and Optimal Control. New Jersey: Prentice Hall, 
1996. 


10 




Figure 1: NASA Flying Test Platform 
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Figure 2: Measured fan speed step responses and predicted fan speed responses (dashed line) by 
the mathematical model. 



Figure 3: Control synthesis interconnection 
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Figure 4: Weighting functions on fan thrust, vane angle, and gimbal torque commands. 






Figure 5: Disturbance forces and moments due to the wind gust. 
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Figure 6: Euler angle responses with the nominal VSVMG wheel speed (1000 rpm). 
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Figure 7: Euler angle responses with different initial VSCMG wheel speeds. 
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Figure 8: Incremental fan speed responses. 






Figure 9: Vane deflection angle responses 
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Figure 10: Girnbal angle responses. 
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Figure 11: VSCMG wheel speed responses at Cl 1000 rpm. 
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Figure 12: VSCMG wheel speed responses at Cl = 2500, 6000, 9000 rpm. 






Figure 13: RMS of gimbal angle, wheel speed, fan speed and vane deflection angle at 1000 rpm 
initial wheel speed. 
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Figure 14: RMS of gimbal angle, wheel speed, fan speed and vane deflection angle. 
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